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ABSTRACT 


This  thesis  examines  the  possibility  of  integrating  an  acoustic  modem  onto  a  Wave  Glider 
with  the  goal  of  providing  a  bidirectional  communications  gateway  for  submerged 
sensors,  platforms,  and  networks.  The  Wave  Glider  unmanned  surface  vehicle 
continuously  harvests  energy  from  the  environment  and  is  able  to  hold  station  without 
needing  to  refuel.  A  unique  two-body  architecture  and  wing  system  directly  converts 
wave  motion  into  thrust,  and  solar  panels  provide  electricity  for  sensor  payloads.  Data 
messages  are  transmitted  to  shore  via  satellite,  and  the  continuous  surface  presence 
means  that  data  can  be  delivered  in  real  time  as  it  is  collected.  The  objective  of  this  thesis 
is  to  identify  the  best  location  for  an  acoustic  modem  on  the  Wave  Glider,  considering  the 
factors  of  hydrodynamic  drag  on  the  vehicle  and  acoustic  performance  of  the  modem. 
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I.  INTRODUCTION 


A,  PROBLEM  STATEMENT 

Emerging  underwater  surveillanee  sensors  and  sensor  networks  will  faeilitate  the 
detection  and  deterrence  of  potential  threats  to  homeland  security.  The  NPS  Seaweb 
project  currently  uses  a  system  of  acoustic  modems  capable  of  underwater  networked 
communication.  A  critical  node  in  the  Seaweb  network  is  the  gateway  node  located  on 
the  ocean  surface  that  connects  the  underwater  nodes  to  command  centers.  Gateway 
nodes  usually  take  the  form  of  a  gateway  buoy.  However,  moored  gateway  buoys  can  be 
expensive  to  maintain  and  are  vulnerable  to  natural  and  human  threats.  This  thesis 
examines  an  alternate  gateway  node  concept,  the  incorporation  of  a  networked  acoustic 
modem  onto  a  versatile  autonomous  marine  craft,  the  Liquid  Robotics  Wave  Glider. 

B,  OBJECTIVE 

The  objective  of  this  thesis  is  to  optimize  the  integration  of  the  acoustic  modem 
onto  the  Wave  Glider.  An  analysis  is  presented  of  the  possible  drawbacks  or  positive 
qualities  associated  with  candidate  locations  for  the  modem.  Considerations  made  in  this 
analysis  include  drag  and  hydrodynamic  performance,  effectiveness  of  acoustic 
communication,  noise  limitations,  shadowing/interference  of  the  vessel,  and  availability 
of  power. 
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II.  BACKGROUND 


A.  SEAWEB  UNDERWATER  ACOUSTIC  NETWORK 

Seaweb  is  a  system  of  acoustic  modems  capable  of  underwater  networked 
communications.  The  network  consists  of  distributed  fixed  and  mobile  nodes  as 
illustrated  in  Figure  1.  The  area  coverage  is  scaleable  according  to  the  environmental 
acoustic  conditions  and  the  number  of  available  nodes.  Typically,  the  network  nodes  are 
expendable,  battery-limited,  and  placed  relatively  close  to  the  sea  floor.  This  is  effective 
in  areas  where  cabled  or  buoyed  sensor  arrays  could  be  subject  to  dangers  such  as 
trawling,  tampering,  and  collisions  with  ship  traffic.  Although,  the  concept  of  Seaweb  has 
been  proven  to  work,  there  remain  aspects  of  the  capability  that  are  undergoing  further 
improvement.  These  include  situation-adaptive,  self-configuring,  and  self-healing 
mechanisms  to  increase  the  longevity  of  Seaweb  components  in  situations  where 
autonomy  is  demanded.  The  aim  of  this  thesis  is  to  improve  the  reliability  and 
independence  of  gateway  nodes  at  the  sea  surface  [1]. 

Information  is  transmitted  between  nodes  acoustically  underwater  until  it  reaches 
a  gateway  node.  The  gateway  node  is  equipped  with  a  modem  that  converts  the  acoustic 
information  into  a  signal  that  is  transmitted  to  an  Iridium  satellite  and  then  redirected  to 
command  centers  monitoring  the  observed  area.  Currently,  gateway  nodes  are  buoys  that 
maintain  station  by  using  an  anchor  and  swivel  system  that  is  moored  to  the  ocean  floor. 
These  gateway  nodes  are  the  critical  interface  between  the  underwater  network  domain 
and  the  above-water  satellite  and  terrestrial  domains.  Ideally,  they  are  located  within  the 
Seaweb  network  so  as  to  maximize  the  number  of  underwater  nodes  that  they  can 
communicate  with  while  also  avoiding  shore  hazards  and  shipping  lanes  [1]. 
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Figure  1 .  The  deseribed  Seaweb  arehiteeture.  Information  is  eolleeted  by  fixed  and 

mobile  nodes,  sent  to  a  gateway  node,  and  then  transmitted  to  a  command 
center  via  the  Iridium  satellite  constellation. 


Currently,  these  buoys  are  statically  moored  and  require  maintenance.  However, 
what  if  the  modems  were  attached  to  an  unmanned  surface  vehicle  (USV)  that  could 
harness  power  from  the  ocean  itself  to  maximize  station  time  and  reduce  network 
maintenance?  Seaweb  engineers  have  been  seeking  an  appropriate  USV  for  this  purpose 
since  2005  [2].  Incorporating  an  acoustic  modem  onto  a  suitable  USV  would  greatly 
extend  the  capabilities  of  an  information  gathering  underwater  network  such  as  Seaweb 
both  in  terms  of  reliability  and  minimization  of  maintenance. 

B.  WAVE  GLIDER  UNMANNED  SURFACE  VEHICLE 

Liquid  Robotics’  Wave  Glider  is  a  wave-powered  USV.  By  using  the  vertical 
motion  of  the  ocean  surface,  it  converts  potential  energy  into  forward  propulsion  in  order 
to  make  way  or  maintain  station.  This  is  achieved  by  using  a  two-body  system  [3].  The 
lower  body  of  the  vehicle  stays  completely  submerged  and  consists  of  six  swiveling  fins 
that  generate  horizontal  forward  thrust  when  the  vehicle  is  moved  vertically  up  and  down. 
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The  upper  body  of  the  Wave  Glider  uses  buoyancy  to  remain  on  the  ocean  surface  and 
pull  the  lower  portion  up  along  with  it  upon  encountering  a  wave.  Using  this  method  of 
propulsion,  the  vehicle  can  achieve  up  to  2  knots  through  the  water.  Clear  advantages  of 
this  energy-harvesting  system  include  lack  of  any  fuel  consumption,  long  endurance,  and 
relatively  little  noise  to  interfere  with  acoustic  communications.  There  is  also  an 
increased  stealth  capability  since  the  main  source  of  propulsion  is  7  meters  below  the 
surface  of  the  water.  The  height  of  the  vehicle  above  the  water  is  no  more  than 
5  to  10  centimeters  [4].  The  manufacturer  of  the  Wave  Glider  refers  to  the  upper  body  as 
the  “float”  and  the  lower  body  as  the  “glider.”  These  terms  will  now  be  used  in  reference 
to  the  two  different  bodies. 


Figure  2.  Potential  energy  from  changes  in  elevation  of  the  Wave  Glider  is  converted 

into  thrust  at  the  glider  portion  of  the  vehicle. 

The  vehicle  is  equipped  with  solar  panels  on  the  float  that  can  provide  up  to  SOW 

of  hotel  power  to  onboard  electronics  or  other  instruments  requiring  power.  Excess  solar 

energy  is  stored  in  a  battery,  as  a  reserve  for  night  operations  and  occurrences  of  high 

demand.  The  solar  energy  system  has  ample  capacity  to  power  an  acoustic  modem  [4]. 
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Figure  3.  A  Wave  Glider  in  operation.  The  float  is  attached  to  the  glider 

via  a  7-meter  tether. 


In  terms  of  survivability  in  harsh  weather  conditions,  the  Wave  Glider  also  excels. 
The  vehicle  was  tested  in  Hurricane  Flossie  in  2007  and  was  able  to  survive  10-foot  seas 
and  40-knot  winds.  While  there  is  a  limit  to  making  way  in  increasingly  rough  sea  states 
because  of  strong  currents  and  multidirectional  waves,  fairly  high  sea  states  will  actually 
help  in  terms  of  maintaining  station,  since  vertical  motion  is  the  vehicle’s  main  source  of 
propulsion  [4]. 

Integrating  an  acoustic  modem  onto  the  Wave  Glider  presents  several  challenges. 
These  include  minimization  of  noise,  reduction  of  drag,  and  identifying  a  preferable 
location  for  acoustic  communications  performance.  The  acoustic  modem  could  be  placed 
underneath  the  upper  body  of  the  Wave  Glider,  on  the  lower  body  of  the  Wave  Glider,  or 
it  could  be  placed  in  a  separate  tow  body,  which  would  be  attached  by  a  tow  to  the  glider. 
Each  of  these  possibilities  has  advantages  and  disadvantages  in  the  considerations  for 
placement  mentioned  earlier,  and  the  goal  was  to  identify  which  location  was  the  most 
advantageous. 
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III.  INTEGRATION  OF  THE  ACOUSTIC  MODEM 


A,  PLACEMENT  OPTIONS 

1,  Mounted  underneath  the  Float 

The  float  portion  of  the  vehiele  is  made  of  fiberglass  and  is  easily  modified.  For 
mounting  an  instrument  on  the  underside  of  the  float,  a  hole  could  be  cut  into  the  body, 
and  the  equipment  could  be  fitted  inside.  This  has  an  advantage  of  minimizing  the  cross- 
sectional  area  of  the  sensor  or  instrument  that  would  be  exposed  to  drag  forces.  However, 
if  the  sensor  is  acoustical  in  nature,  as  is  our  modem,  it  has  the  disadvantage  of  being 
close  to  the  surface.  This  makes  it  more  susceptible  to  noise.  Another  advantage  is  the 
direct  access  to  power  and  to  the  onboard  electronics  in  the  payload  bays  within  the  float. 

Some  instruments  have  previously  been  integrated  into  the  Wave  Glider  float.  In 
one  experiment.  Teledyne  Benthos  and  Liquid  Robotics  were  interested  in  integrating  an 
Acoustic  Doppler  Current  Profiler  (ADCP)  with  the  vehicle  in  order  to  collect  current 
measurements  near  the  surface  of  the  water.  A  hole  was  cut  in  the  hull  of  the  float  and  the 
ADCP  was  mounted  in  it.  A  fairing  was  also  constructed  around  the  ADCP  in  order  to 
reduce  the  drag  associated  with  it.  This  integration  can  be  seen  in  Figure  4  [5]. 


Figure  4.  An  Acoustic  Doppler  Current  Profiler  mounted  to  the  float. 

The  streamlined  shape  was  used  to  reduce  drag. 
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2. 


Mounted  on  the  Glider 


Mounting  the  modem  on  the  glider  has  the  advantage  of  being  away  from  ocean 
turbulence,  bubbles,  and  noise.  However,  unlike  the  float,  there  is  no  existing  volume  for 
encapsulating  the  modem  and  hence  reducing  the  cross-sectional  area  that  would  be 
presented  to  the  flow  field.  Thus,  the  amount  of  drag  created  due  to  the  modem  would 
increase  in  this  scheme.  Also,  all  of  the  vehicle’s  propulsion  is  generated  from  the  glider 
portion  of  the  vehicle.  Attaching  more  mass  onto  the  glider  without  accounting  for  the 
shift  in  the  center  of  gravity  and  other  hydrodynamic  properties  of  the  vehicle  might 
severely  impact  the  efficiency  of  the  propulsion  system.  Lastly,  wiring  would  have  to  be 
implemented  in  order  to  connect  the  modem  to  the  instrumentation  and  power  source 
located  on  the  float. 

Integration  schemes  where  instrumentation  is  attached  to  the  glider  section  have 
already  been  used  in  multiple  experiments.  While  thus  far  the  integration  has  been  crude 
with  respect  to  minimization  of  drag,  usually  it  meets  the  needs  for  the  experiment  being 
carried  out  at  the  time.  One  such  integration  is  shown  in  Figure  5.  In  this  experiment,  an 
active  transducer  and  a  passive  hydrophone  were  placed  on  the  tail  end  of  the  glider 
body.  No  effort  was  made  to  streamline  the  transducer  or  integrate  it  into  the  structure  of 
the  glider  [6]. 


Figure  5.  An  experiment  in  which  an  active  transducer  and  a  passive  hydrophone 

were  mounted  onto  the  glider. 
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3. 


Pulled  in  a  Tow  Body 


The  tow-body  configuration  has  many  appealing  characteristics  about  it  compared 
with  the  latter  two  configurations.  It  has  both  the  advantages  of  being  away  from  the 
surface  of  the  water,  and  being  completely  customizable  in  terms  of  shape  in  order  to 
minimize  drag.  One  would  think  that  the  extra  cable  connecting  the  tow  body  to  the 
glider  would  produce  a  significant  drag  on  its  own.  However,  a  mechanically  compliant, 
low-drag  cable  has  been  developed  to  mitigate  this  issue.  The  biggest  disadvantage  with  a 
tow  body  is  the  increased  handling  complexity.  Deployment  and  recovery  operations  at 
sea  will  be  more  difficult  and  more  risky  with  a  three-body  configuration. 

In  March  2012,  Cornell,  BioSonics,  and  Liquid  Robotics  worked  together  to 
perform  some  tests  on  a  BioSonics  dual-frequency  echosounder.  In  this  experiment,  a 
Wave  Glider  pulled  a  tow  body  housing  a  70  kHz  and  200  kHz  transducer.  In  Figure  6, 
we  can  see  the  tow  body  [7]. 


Figure  6.  The  Wave  Glider  and  the  tow  body  used  in  the  dual-frequency  echosounder 

experiment  off  the  coast  of  Hawaii. 
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4,  Separation  of  Receiver  and  Transmitter 

Another  possibility  that  was  considered  but  not  analyzed  in  this  thesis  is  that  of 
separating  the  transmitter  and  receiver  portions  of  the  modem.  The  larger  transmitter 
could  be  placed  underneath  the  float  and  gain  the  benefit  of  minimized  cross-sectional 
area.  The  receiver  would  be  placed  on  the  glider  or  in  a  towed  hydrophone  array  to  avoid 
surface  noise.  This  complicates  the  integration  because  of  the  separation  of  modem  parts, 
additional  wiring  required,  and  delivery  of  power.  However,  it  should  still  be  examined  in 
future  research. 

B,  INTEGRATION  CONSIDERATIONS 

1.  Drag 

The  complete  configurability  of  the  tow  body  means  that  we  could  optimize  the 
shape  to  reduce  its  associated  drag.  It  will  be  informative  to  estimate  how  much  drag 
would  be  produced  from  placing  the  transducer  on  either  the  glider  body  or  the  float  in 
order  to  quantify  the  advantage  of  placing  the  modem  in  the  tow  body.  In  Chapter  IV,  a 
description  is  given  of  the  approximations  and  theory  used  to  provide  an  estimate  of  this 
advantage. 

2.  Acoustic  Communications 

It  would  also  be  helpful  to  identify  what  impact  the  depth  of  the  sensor  would 
have  on  the  acoustic  propagation  path  to  a  node  below  it.  Also  of  interest  is  whether  or 
not  a  modem  mounted  on  the  float  would  encounter  problems  due  to  the  glider  being 
beneath  the  sensor  at  some  angle.  In  Chapter  V,  an  analysis  of  ray  theory  is  given  along 
with  some  propagation  modeling  using  Bellhop  and  Matlab. 

3.  Noise 

Lastly,  we  wish  to  estimate  how  much  noise  is  associated  with  each  modem 
location  at  the  frequencies  of  interest.  As  mentioned  previously,  the  float  configuration 
would  be  subject  to  surface  noise  that  would  be  lessened  by  placing  the  modem  deeper  in 
the  water.  An  analysis  of  different  sources  of  noise  and  their  impact  on  the  system  is 
provided  in  Chapter  VI. 
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IV.  DRAG  CONSIDERATIONS 


A.  BACKGROUND 

1.  Drag  Coefficient 

An  object  moving  through  a  fluid  at  any  velocity  experiences  a  resistive  force 
known  as  drag.  Drag  is  caused  by  the  pressure  and  shear  stress  distributions  on  an  object 
due  to  a  particular  fluid’s  viscosity  and  flow  characteristics.  If  the  pressure  and  shear 
stress  distributions  over  the  object  are  known,  one  would  simply  integrate  the  pressure 
and  shear  stress  distributions  over  the  surface  area  of  the  object  in  order  to  obtain  the  net 
force  acting  in  the  opposite  direction  of  the  velocity.  However,  most  of  the  time,  these 
distributions  are  not  known,  and  empirical  data  must  be  used  in  order  to  solve  for  the 
drag  forces  on  the  object.  Typically,  this  is  done  by  using  a  dimensionless  quantity 
known  as  the  drag  coefficient.  The  drag  coefficient  is  defined  by 


(4.1) 


Fq 

\pAV^ 


Here  Fd  is  equal  to  the  drag  force,  p  is  the  density  of  the  fluid,  V  is  the  velocity  of 
the  object  relative  to  the  fluid,  and  A  is  the  characteristic  area  of  the  object  (usually  the 
frontal  area  or  the  surface  normal  to  the  direction  of  flow)  [8]. 

Sometimes  the  drag  coefficient  is  separated  into  two  component  coefficients. 

(4.2) 

The  first  component  is  the  coefficient  of  drag  due  to  friction  and  the  second  is  the 
coefficient  of  drag  due  to  pressure.  Friction  effects  are  more  noticeable  when  the 
viscosity  of  the  fluid  is  high  and  there  is  a  lot  of  shear  stress  acting  on  the  object  while 
the  pressure  drag  is  more  strongly  dependent  on  the  shape  or  form  of  the  object. 
Depending  on  the  fluid  that  an  object  is  in,  the  drag  coefficient  can  have  very  different 
values  [9]. 
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To  assess  the  increased  drag  caused  by  adding  an  acoustic  modem  onto  the  Wave 
Glider,  we  approximate  the  acoustic  modem  as  a  cylinder  in  incompressible  fluid  flow.  It 
should  be  noted  that  this  analysis  applies  to  the  case  where  the  modem  is  mounted 
without  any  fairing. 


2.  Drag  Coefficient  of  a  Cylinder  in  Viscous  Incompressible  Flow 

Most  information  relating  to  drag  has  been  attained  by  conducting  numerous  flow 
experiments  on  the  object  of  interest.  The  exact  size  and  shape  of  the  modem  we  plan  to 
use  are  not  available.  However,  we  can  approximate  its  dimensions  and  draw  some 
conclusions  as  to  its  drag  characteristics. 

We  approximate  the  acoustic  modem  as  a  cylinder.  The  drag  coefficient  for  a 
cylinder  in  incompressible  flow  normal  to  the  axis  can  be  found  in  a  table  or  book.  The 
drag  coefficient  for  a  cylinder  in  incompressible  flow  is 


(4.3) 


■  ^Df  ^Dp 


5.93 


+  1.17 


Here,  the  drag  coefficient  is  expressed  as  a  function  of  the  Reynolds  number.  The 
Reynolds  number  is  a  dimensionless  quantity  that  describes  the  ratio  of  inertial  effects  to 
viscous  effects  [9]. 


(4.4) 


Re  = 


pVL 


Here,  p  is  the  density  of  the  fluid,  V  is  the  velocity  of  the  object  relative  to  the 
fluid,  L  is  the  traveled  length  of  the  fluid,  and  p  is  the  dynamic  viscosity  of  the  fluid.  So 
in  order  to  figure  out  the  coefficient  of  drag,  we  must  first  calculate  the  Reynolds  number 
of  seawater.  Plugging  in  all  of  the  values  for  seawater,  assuming  that  the  traveled  length 
of  the  fluid  is  half  the  circumference  of  our  cylinder,  and  assuming  a  velocity  of  1  knot 
we  get  a  Reynolds  number  of  77000  (a  large  Reynolds  number).  This  is  to  be  expected, 
because  seawater  is  not  very  viscous.  Therefore,  there  is  negligible  friction  drag  in  our 
analysis.  Pressure  drag  dominates,  and  we  are  able  to  make  assumptions  based  on  very 
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high  Reynolds  number  flows.  If  we  look  at  our  coeffieient  of  drag  now,  the  frietion 
component  is  negligible,  because  the  Reynolds  number  is  so  large. 

(4.5)  c^=1.17 

This  also  assumes  though  that  the  cylinder  is  infinitely  long  in  a  2  dimensional 
fluid  flow.  In  Figure  7,  we  can  see  that  the  coefficient  of  drag  of  a  cylinder  depends  on 
the  aspect  ratio  of  the  cylinder  [10].  For  our  modem,  let  us  assume  that  the  length  of  the 
cylinder  is  approximately  1.5  times  that  of  the  diameter. 


Figure  7.  A  graph  relating  the  aspect  ratio  of  a  cylinder  with  two  free  ends  to  its  drag 
coefficient.  We  can  see  that  as  the  length  of  the  cylinder  becomes  greater 
than  its  diameter,  the  drag  coefficient  approaches  the  value  1.17. 


Thus  our  new  coefficient  of  drag  becomes  smaller. 
(4.6)  c^=0.7 
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3,  Drag  Force  of  Acoustic  Modem  in  Seawater 

Now  that  we  have  an  approximate  eoeffieient  of  drag  for  the  modem  we  ean  get 
an  approximation  for  the  drag  foree  exerted  on  it  while  it  moves  through  the  water.  The 
previously  mentioned  equation  of  drag  ean  be  rearranged  to  solve  for  the  foree. 

1  7 

(4.7)  F^=-pc^AV 

If  we  plug  in  values  on  the  right  side  of  the  equation,  we  obtain  the  drag  foree 
assoeiated  with  adding  a  modem  to  the  vehiele,  depending  on  the  amount  of  exposed  area 
and  the  speed  at  whieh  the  vehiele  is  traveling.  However,  this  does  not  help  us  ealeulate 
the  reduetion  in  speed  experieneed  by  the  vehiele  until  we  know  the  amount  of  thrust 
generated  by  the  propulsive  system  of  the  Wave  Glider. 

B,  ANALYSIS  OF  WAVE  GLIDER  PROPULSION  SYSTEM 

1.  Previous  Work 

Work  has  been  done  previously  quantifying  the  foree  of  thrust  of  a  wave 
propelled  vehicle  as  a  function  of  wave  height  and  period  [11].  Using  a  wave  tank  and  a 
cable  attached  to  a  force  meter,  one  can  measure  the  tension  in  the  cable  and  obtain 
values  for  the  thrust  of  the  vehicle  as  the  wave  height  and  period  in  the  tank  is  adjusted. 
This  approach  was  considered  in  this  analysis  as  viable  method  of  determining  a 
reduction  in  speed  of  the  Wave  Glider.  However,  here  we  limit  ourselves  to  a  theoretical 
estimate  based  on  the  power  of  ocean  waves  and  power/thrust  relations. 

2,  Power  Generated  by  Ocean  Waves 

The  amount  of  thrust  generated  by  the  Wave  Glider  is  dependent  on  the  changes 
in  potential  energy  associated  with  the  vertical  displacement  of  the  float.  This  is 
converted  into  forward  thrust  by  the  flapping  motion  of  the  fins  on  the  glider  body.  The 
theoretical  amount  of  power  in  a  monochromatic  deep  ocean  wave  is  a  function  of  the 
wave  height  H,  and  the  wave  period  T. 

(4.8) 

14 


Here  p  is  the  density  of  seawater,  and  g  is  the  aeeeleration  due  to  gravity.  This  is 
aetually  the  power  of  the  wave  per  meter  of  erest  length  [12],  Crest  length  is  a 
measurable  parameter  related  to  the  speetral  density  of  the  oeean  waves  [13].  Many 
studies  have  gone  into  produeing  distributions  of  wave  erest  length  in  the  oeean  under 
different  weather  and  sea  eonditions.  For  the  purposes  of  this  study,  I  will  approximate 
the  length  of  the  wave  erest  as  the  speed  at  which  the  Wave  Glider  is  traveling  multiplied 
by  the  wave  period,  since  the  order  of  wave  crests  length  is  comparable  to  that  of  the 
wave  length  [14].  We  now  have  a  new  equation  for  the  power  of  each  wave. 


(4.9) 


n 


wave 


pg^H^T^V 

^In 


[W] 


The  vehicle  is  of  course  not  able  to  use  all  of  this  power  for  thrust  as  nothing  is 
100%  efficient.  The  next  step  in  finding  the  thrust  generated  is  to  estimate  the  efficiency 
of  the  vehicle. 


3,  Wave  Glider  Propulsion  Efficiency 

Flapping  foil  motion  used  for  propulsion  is  actually  very  efficient  when  compared 
to  other  methods  of  propulsion.  Much  research  is  still  being  conducted  on  how  certain 
birds,  fish,  and  insects  are  able  to  swim/fiy  so  much  more  efficiently  than  some  of  the 
propellers  in  use  today.  The  reason  it  is  not  used  for  everything  is  because  it  is  somewhat 
limited  in  the  amount  of  thrust  that  it  can  produce.  Higher  levels  of  thrust  can  be 
generated;  however,  much  of  the  efficiency  is  sacrificed.  The  efficiency  of  the  flapping 
foil  can  be  characterized  as  a  function  of  the  Strouhal  number,  the  heave  amplitude,  the 
angle  of  attack  of  the  foil,  and  the  phase  angle  between  the  pitch  and  the  heave.  In 
Figure  8,  a  graph  shows  the  different  contours  of  the  efficiency  of  a  flapping  foil  as  a 
function  of  the  angle  of  attack  and  the  Strouhal  number  [15]. 
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Figure  8.  A  contour  plot  of  the  efficiency  of  an  airfoil  with  a  given  ratio  between  the 

heave  amplitude  and  the  chord  length  c 


Using  some  of  the  empirical  data  from  this  particular  study,  we  can  find  an 
approximation  to  the  propulsive  efficiency  of  the  Wave  Glider  once  we  know  values  for  a 
couple  of  the  parameters.  The  first  one  we  can  find  is  the  Strouhal  number.  The  Strouhal 

number  is  a  function  of  the  heave  amplitude  ,  the  angular  frequency  of  the  flapping  foil 
CO ,  and  the  velocity  of  the  foil  V  (or  in  this  case  the  speed  of  the  Wave  Glider)  [14]. 


(4.10) 


_  ^crh^co 
V 


We  can  calculate  the  Strouhal  numbers  for  our  vehicle  by  using  average  speeds  in 
various  sea  states.  These  sea  states  also  correspond  to  the  wave  heights  experienced  by 
the  vehicle  and  therefore  to  the  heave  amplitude.  Lastly,  the  angular  frequency  with 
which  the  foils  oscillate  should  be  related  to  the  period  of  the  waves  according  to 
(x)=1kIT.  Typically,  wave  periods  in  the  ocean  range  from  6-20  seconds  [12],  so  let  us  use 
a  value  of  T  =  6  seconds  for  the  wave  period  and  thus  the  angular  frequency  of  the 
oscillating  foils  is  7i/3  [rad/s].  Table  1  shows  the  average  speed  of  the  wave  glider 
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in  varying  sea  states  [4],  I  have  adjusted  the  table  so  that  it  also  shows  the  wave  heights 
corresponding  to  those  sea  states  (using  the  Douglas  Sea  Scale)  and  the  average  speed  in 
meters  per  second. 


Sea  State 

Wave  Height 
(m) 

Wave  Glider 
Speed  (knots) 

Wave  Glider 
Speed  (m/s) 

Oscillating 

Frequency 

(rad/s) 

Strouhal 

Number 

Flat  Calm 

0.00 

0.00 

0.00 

1.05 

0.00 

Sea  State  0 

0.01 

0.38 

0.19 

1.05 

0.68 

Sea  State  1 

0.05 

1.00 

0.51 

1.05 

1.28 

Sea  State  2 

0.30 

1.63 

0.84 

1.05 

4.72 

Sea  State  3 

0.88 

1.88 

0.96 

1.05 

11.94 

Table  1 .  Data  relating  the  speed  of  the  Wave  Glider  to  the  sea  state.  A  wave  period 
of  6s  is  assumed  thus  producing  an  oscillating  frequency  of  ®  =  1.05  rad/s 
Strouhal  numbers  were  calculated  from  each  corresponding  heave 
amplitude  (wave  height),  speed  of  the  wave  glider,  and  oscillating 

frequency. 

Unfortunately,  the  study  of  efficiency  of  flapping  foils  did  not  go  to  higher 
Strouhal  numbers  [15],  and  the  ratio  between  the  heave  amplitude  and  the  chord  length 
was  kept  to  1 .  Although  our  heave  amplitudes  are  much  greater  than  the  chord  length  of 
our  foils,  I  decided  to  extrapolate  from  their  data  to  estimate  the  efficiency  of  the  vehicle. 
The  Strouhal  number  derived  from  the  sea  state  0  value  in  Table  1  is  0.68  and  the  angle 
of  attack  of  the  foils  of  the  wave  glider  is  around  40°.  Using  these  values  and  the  contour 
plot  in  Figure  8,  I  approximated  the  net  efficiency  of  the  wave  glider  to  be  around 
77  =  0.45 . 

4,  Wave  Power  Utilized  by  the  Wave  Glider 

Now  that  we  have  a  rough  estimate  of  the  efficiency  of  the  Wave  Glider  and  the 

power  contained  within  an  ocean  wave  of  a  given  height,  we  can  approximate  the  power 

available  for  forward  thrust.  We  simply  multiply  the  efficiency  of  the  system  by  the 

power  of  the  wave  to  get  the  motive  power  of  the  Wave  Glider.  Table  2  takes  the  values 

for  wave  height  from  Table  1  and  lists  the  power  associated  with  these  wave  heights  by 
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using  Equation  (4.8)  (once  again  assuming  a  wave  period  of  6  seeonds).  These  values  are 
then  multiplied  by  the  effieieney  of  the  Wave  Glider  to  get  the  power  available  for  thrust. 


Sea  State 

Wave  Height 
(m) 

Wave  Glider 
Speed  (m/s) 

Power  of 
Ocean  Waves 
(W) 

Efficiency 
of  Wave 
Glider  rj 

Power 
Available  for 
Thrust  (W) 

Flat  Calm 

0.00 

0.00 

0.00 

0.45 

0.00 

Sea  State  0 

0.01 

0.19 

0.680 

0.45 

0.306 

Sea  State  1 

0.05 

0.51 

45.3 

0.45 

20.4 

Sea  State  2 

0.30 

0.84 

2650 

0.45 

1190 

Sea  State  3 

0.88 

0.96 

26000 

0.45 

11700 

Table  2.  Values  for  the  Wave  Glider’s  available  thrust  are  ealeulated  using  the  power 
available  in  an  ocean  wave  of  the  indicated  height  and  the  approximate 
effieieney  of  the  Wave  Glider. 


5.  Drag  and  Thrust  Forces  of  the  Wave  Glider 

Using  some  basie  meehanies,  we  ean  now  ealeulate  the  drag  and  thrust  foree  of 
the  Wave  Glider  by  using  the  power  values  that  we  just  ealeulated.  For  propulsion  there 
is  a  relation  between  the  thrust  force,  veloeity  of  the  vehicle,  and  the  power  of  the 
propulsion  [16]. 

(4-11)  ^WG~^Thrusy 

Assuming  that  the  veloeity  of  the  Wave  Glider  is  approximately  eonstant  in  these 
sea  state  conditions,  we  divide  the  power  by  the  speed  of  the  Wave  Glider  to  find  the 
thrust  foree.  We  should  also  note  though  that  sinee  the  vehiele  is  going  at  a  eonstant 
speed  and  not  aoeelerating,  there  ean  be  no  net  forees  aeting  on  the  vehiele.  Figure  9 
shows  a  free-body  diagram  of  the  Wave  Glider  [6]. 
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This  means  that  the  magnitude  of  the  force  of  drag  must  be  equal  to  the  thrust 
force.  This  is  important,  because  based  on  the  speed  it  is  going  we  can  use  Equation  (4.7) 
to  calculate  the  coefficient  of  drag. 


(4.12) 

o 

II 

II 

(4.13) 

If  I  = 

1  Thrust  \ 

F, 


Drag 


6.  Coefficient  of  Drag  of  the  Wave  Glider 

Since  we  know  the  thrust  force,  we  can  set  it  equal  to  the  drag  force.  We  then  use 
Equation  (4.1)  to  solve  for  the  drag  coefficient.  Table  3  shows  the  values  for  the  thrust 
and  drag  forces  and  the  coefficient  of  drag.  The  frontal  area  of  the  vehicle  was 
approximated  at  0.060  by  using  dimensions  provided  by  a  Eiquid  Robotics 
specifications  sheet  and  some  measurements  of  my  own.  These  area  measurements 
assumed  an  ideal  case  in  which  the  Wave  Glider  was  perfectly  level. 
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Sea  State 

Power 

Available  for 
Thrust  (W) 

Wave  Glider 
Speed  (m/s) 

Thrust  and 
Drag  Force(N) 

Coefficient  of 
Drag  of  Wave 
Glider 

Flat  Calm 

0.00 

0.00 

0.00 

0.00 

n/a 

Sea  State  0 

0.306 

0.19 

1.59 

1.14 

1.39 

Sea  State  1 

20.4 

0.51 

39.7 

8.14 

4.87 

Sea  State  2 

1190 

0.84 

1430 

21.5 

66.4 

Sea  State  3 

11700 

0.96 

12100 

28.6 

425 

Table  3.  Estimates  for  the  Coeffieient  of  Drag  of  the  Wave  Glider 


The  eoeffieient  of  drag  should  be  relatively  eonstant  sinee  the  only  parameters 
that  are  ehanging  in  our  model  are  the  power  of  the  wave  and  the  speed  of  the  wave 
glider.  Unfortunately,  as  we  ean  see  from  Table  3,  our  results  are  not  eonsistent  with  this 
expeetation.  The  disagreement  probably  results  from  us  using  an  overly  simplistie  model 
for  the  power  of  eaeh  wave  and  assuming  that  the  effieieney  of  the  vehiele  is  eonstant  as 
the  wave  height  and  the  speed  of  the  Wave  Glider  inerease.  Also,  wave  heights  assoeiated 
with  different  sea  states  ean  vary  quite  a  bit.  Although  some  of  our  estimates  for  the 
eoeffieient  of  drag  seem  quite  far  off,  our  estimates  for  the  eoeffieient  of  drag  resulting 
from  the  data  assoeiated  with  sea  state  0  seem  to  be  fairly  reasonable.  Thus  we  will  use 
the  eoeffieient  of  drag  value  from  sea  state  0  in  all  of  our  future  ealeulations. 

(4.14)  c^=1.39 

C.  CALCULATING  REDUCTION  IN  WAVE  GLIDER  SPEED 

1.  Using  Power,  Drag,  and  Thrust  Relations  to  Find  Speed  Reduction 

To  find  the  reduetion  in  speed,  we  start  with  Equation  (4.11)  and  realize  onee 
again  that  the  thrust  foree  will  be  equal  to  the  drag  foree  onee  the  Wave  Glider  reaehes  a 
eonstant  speed. 

(4.15)  = 
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This  drag  force  is  the  sum  of  drag  forces  from  the  modem  or  the  tow  body,  the 
drag  force  from  the  umbilical  cable  and  the  drag  forees  on  the  Wave  Glider. 

(4-16)  )  V 

We  then  also  plug  Equation  (4.7)  into  Equation  (4.16)  for  eaeh  of  the  drag  forces. 


(4.17) 


^IVG 


1 


1 


1 


A 


P^WG^  WG^  ^  P^Modem^ModettX  ^  P^Cable^Cable^ 

1  1  J 

1 


V 


(4.18) 


We  can  factor  out  —  py^ . 

2 


1 


^WG  P^  ipWG^WG  *^Modem^Modem'^  ^Cable^Cahle) 


Now  we  can  solve  for  the  new  velocity  of  the  vehicle  with  the  modem  added, 
provided  that  we  also  know  the  coeffieients  of  drag  and  eross-sectional  areas  of  both  the 
modem  and  the  Wave  Glider. 


(4.19) 


V  = 


2n 


n(i/3) 


WG 


Pi 


^WG^WG  ^Modern  ^Modern  ^ Cable  ^Cable  y 


2,  Wave  Glider  Velocity  after  Mounting  underneath  the  Float 

If  the  modem  is  mounted  underneath  the  float  without  any  fairing,  most  of  the 
cross-sectional  area  is  not  exposed  to  the  flow  field.  Since  we  assumed  earlier  that  the 
modem  is  about  1.5  times  as  long  as  its  diameter,  let  us  use  a  length  of  15  cm  and  a 
diameter  of  10  cm  to  approximate  the  modem.  Sinee  we  are  analyzing  the  situation  where 
it  is  mounted  on  the  float,  let  us  also  assume  that  ki  of  the  modem  is  exposed.  This  makes 
the  frontal  area  of  the  modem  equal  to  a  plane  of  3.75  cm  by  10  cm.  Eigure  10  shows  the 
geometry. 
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15cm 


10cm 


3.75cm 

i 


10cm 


10cm 


Figure  10.  Twenty- five  pereent  of  the  Cylinder  is  exposed.  Also  the  frontal  area  of  the 
eylinder  ean  be  eonsidered  to  be  a  plane  of  equal  dimensions. 


Now  we  have  frontal  areas  and  eoeffieients  of  drag  for  both  the  wave  glider  and 
the  partially  exposed  modem.  Also,  the  drag  eoeffieient  of  the  faired  eable  is  0.12  [9]  and 
the  eross-seetional  area  is  0.07m  .  Table  4  shows  these  values  and  the  resulting  value  of 

{^WG^WG  ^ Modem^Modem  ^ Cable^Cahle)  ' 


{^WG-^WG  ^Cable-^Cahle) 

Area  of 
Modem  (m^) 

Cj^  of  Modem 

(<^WG^WG  ^Modern ^Modern  ^ Cab/e ^Cab/e) 

0.0916 

0.00375 

0.7 

0.0942 

Table  4.  New  values  for  the  frontal  surfaee  area  and  drag  eoeffieients  are  ealculated 

for  a  float-mounted  modem. 


To  find  the  reduction  in  speed,  we  use  Equation  (4.19),  the  values  for  power 
available  for  thrust  in  Table  3,  and  the  final  value  for  J 

in  Table  4.  As  mentioned  previously,  sinee  we  are  using  the  eoeffieient  of  drag  value  for 
the  Wave  Glider  eorresponding  to  sea  state  0,  we  should  also  use  the  power  and  initial 
speed  values  for  sea  state  0  as  well.  The  3’^‘^  eolumn  of  Table  5  is  the  result  of  using 
Equation  (4.19)  to  find  the  new  resultant  speed,  and  the  4th  eolumn  indieates  the 
percentage  reduction  of  speed. 
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Power 

Available  for 
Thrust  (W) 

Speed 
Without 
Modem  (m/s) 

Calculated 
Speed  With 
Modem  (m/s) 

%  Reduction 
in  Speed 

0.306 

0.193 

0.185 

4.11% 

Table  5.  The  reduetion  in  speed  of  a  float-mounted  modem  at  sea  state  0. 


A  4.11%  reduetion  in  speed  is  not  a  tremendous  penalty.  This  eould  be  redueed 
further  if  a  fairing  was  construeted  around  the  modem.  If  our  main  eoneern  for 
integration  is  drag,  a  float  mounted  modem  is  a  viable  option. 

3,  Wave  Glider  Velocity  after  Mounting  on  Glider 

To  calculate  the  reduction  in  speed  of  mounting  the  modem  to  the  glider,  we  use 
the  same  process  illustrated  in  the  last  section.  The  only  difference  is  that  now  the  modem 
is  considered  to  be  fully  exposed  to  the  water.  The  frontal  area  of  the  cylinder  is  now  a 
plane  of  15  cm  by  10  cm.  The  geometry  is  shown  by  Figure  11. 


15cm 


10cm 


10cm 


Figure  1 1 .  This  time,  the  modem  is  fully  exposed.  Once  again,  the  frontal  area  of  the 

cylinder  can  be  considered  to  be  a  plane  of  equal  dimensions. 
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New  values  for  the  area  and  coefficient  of  drag  are  reflected  in  Table  6. 


{^WG^WG  ^Cable^Cahle) 

Area  of 
Modem  (m^) 

Cjy  of  Modem 

(^WG^WG  ^Modern ^Modern  ^Cab/e^Cah/e) 

0.0916 

0.015 

0.7 

0.102 

Table  6.  New  areas  and  coefficients  of  drag  are  calculated  for  the  glider-mounted 

modem. 


Finally,  the  value  of  +  ^cabie^chie)  is  used  to  calculate  the 

reduction  in  speed.  Once  again  I  used  the  power  values  and  initial  speed  values  from  sea 
state  0.  Table  7  shows  the  calculated  reduction  in  speed  from  a  glider  mounted  modem. 


Power 

Available  for 
Thrust  fW) 

Speed 
Without 
Modem  fm/s) 

Calculated 
Speed  With 
Modem  fm/s) 

%  Reduction 
in  Speed 

0.306 

0.193 

0.180 

6.64% 

Table  7.  The  reduction  in  speed  of  a  glider-mounted  modem  at  sea  state  0. 


While  still  not  a  tremendous  reduction  in  speed,  it  is  still  more  than  the  float 
mounted  case.  However,  it  should  be  noted  that  this  does  not  take  into  account  changes  in 
the  efficiency  of  the  vehicle  due  to  the  modem  being  mounted  on  the  same  portion  of  the 
vehicle  that  provides  all  of  the  propulsion.  More  serious  reductions  in  speed  could 
possibly  occur  if  the  modem  were  not  placed  near  the  center  of  gravity  of  the  vehicle  or  if 
it  somehow  inhibited  the  Wave  Glider’s  method  of  harvesting  wave  energy. 

4,  Wave  Glider  Velocity  after  Adding  a  Tow  Body 

Most  tow  bodies  being  developed  today  have  a  bullet  shaped  profde  in  order  to 
minimize  drag.  Figure  12  shows  an  approximate  shape  for  a  tow  body. 
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Figure  12.  The  side  and  front  profiles  for  an  approximate  tow  body 


We  will  eonsider  a  tow  body  with  a  diameter  that  is  big  enough  to  house  the 
modem  with  D  =18  em.  The  eoeffieient  of  drag  for  sueh  a  shape  is  approximately 
Cjj  =  0.30  and  the  frontal  area  is  the  area  of  a  eirele  of  diameter  D  [9].  Using  the  same 

process  as  the  previous  two  sections,  we  calculate  the  reduction  in  speed.  Table  8  shows 
the  new  total  frontal  areas  and  drag  coefficients. 


{^WG^WG  ^Cable^Cable) 

Area  of  Tow 
Body  (m^) 

Cj^  of  Tow 
Body 

{^WG^WG  ^ Modem^Modem  ^Cable^Cable) 

0.0916 

0.025 

0.3 

0.0991 

Table  8.  Frontal  area  and  coefficient  of  drag  values  after  adding  a  tow  body 


Again,  these  values  are  then  used  to  calculate  a  reduction  in  speed.  Table  9  shows 
the  results. 


Power 

Available  for 
Thrust  fWl 

Speed 
Without 
Modem  fm/s) 

Calculated 
Speed  With 
Modem  fm/s) 

%  Reduction 
in  Speed 

0.306 

0.193 

0.182 

5.7% 

Table  9.  The  reduction  in  speed  of  a  modem  mounted  inside  a  tow  body. 
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This  configuration  lies  in  between  the  float-mounted  and  glider-mounted  cases. 
The  tow  body  may  impaet  the  Wave  Glider’s  propulsive  system.  In  addition,  small 
amounts  of  drag  could  be  induced  by  the  eable  connecting  the  tow  body  to  the  glider,  but 
the  hydrodynamics  of  the  tow  body  can  be  improved  by  incorporating  compliance  and 
fairing  into  the  tow  eable. 

D,  SUMMARY 

From  our  analysis,  the  best  integration  in  terms  of  drag  is  a  modem  mounted 
underneath  the  float  body.  The  reductions  in  speed  for  each  case  are  shown  in  Table  10. 


Integration  Scheme 

%  Reduction  in  Speed 

Float-Mounted 

4.1% 

Glider-Mounted 

6.6% 

Tow-Body  Mounted 

5.7% 

Table  10.  Comparative  reductions  in  speed  for  the  3  different  integration  schemes 

As  mentioned  previously  the  float  mounted  reduction  in  speed  could  be  improved 
by  construeting  a  fairing  around  the  modem.  It  would  seem  that  the  glider  mount  is 
disadvantageous  in  all  respects.  Both  the  float  mount  and  tow-body  mount  avoid  the  issue 
of  interfering  with  the  Wave  Glider’s  propulsion  and  have  favorable  values  for  reduction 
of  speed. 

It  has  been  noted  that  these  values  seem  low  eompared  to  test  results  already 
earried  out.  In  one  particular  test  in  a  Liquid  Robotics  test  facility  in  Kwaihea,  Hawaii, 
the  Wave  Glider  suffered  a  20%  vehicle-speed  reduetion  compared  to  an  unloaded 
vehiele  [17].  Possible  errors  in  our  analysis  may  arise  due  to  our  approximated  oeean 
eonditions,  namely  the  wave  period  and  wave  erest  lengths.  Regardless,  the  hierarehy  of 
the  values  that  we  calculated  seems  reasonable. 
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V.  ACOUSTIC  COMMUNICATIONS 


A.  RAY  PROPAGATION  THEORY 
1.  The  Wave  Equation 

Similarly  to  waves  in  other  media,  acoustic  waves  of  small  enough  amplitude  in 
water  obey  the  linear  wave  equation.  In  terms  of  pressure,  p,  the  spatial  coordinates  x,  y, 
and  z,  the  wave  speed  c,  and  time  t,  the  wave  equation  can  be  expressed  as 


(5.1) 


d^p  _  2  ^ d^p  d^p  d^p^ 
dt^  ^  dx^  dy^  dz^  ^ 


One  way  to  solve  the  wave  equation  is  called  ray  theory.  It  is  useful  for  understanding  an 
acoustic  wave’s  phase  and  “position”  in  terms  of  space  and  time.  Although  waves  cannot 
actually  have  a  position,  it  is  much  easier  to  visualize  how  a  wave  propagates  in  a 
particular  medium  by  using  ray  theory  [18]. 

Without  going  into  too  much  detail,  we  will  briefly  describe  how  one  goes  from 
the  wave  equation  to  ray  theory.  First,  we  use  separation  of  variables  to  go  from  the  wave 
equation  to  the  time-independent  Helmholtz  equation.  Note  that  the  spatial  derivatives  in 
Equation  (5.1)  have  been  condensed  into  and  x,  y,  and  z  into  r. 

(5.2)  p{^r^-0 


Next,  we  plug  in  an  ansatz  solution  to  the  Helmholtz  equation. 


(5.3) 


Here  A{r)  represents  the  amplitude  of  the  wave  and  cox  represents  the  phase  of  the 
wave.  After  some  amount  of  math  and  taking  the  limit  as  O  ^  go  ,  we  arrive  at  a  form  of 
the  eikonal  equation. 

Iy7  P  ^ 

(5.4)  |Vr|  =- 
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From  this  we  can  also  derive  the  transport  equation.  This  describes  how  the 
amplitude  changes  over  the  regime  of  the  phase  front. 

(5.5)  2VA-Vt  +  AV^t  =  0 

The  transport  equation  and  the  eikonal  equation  can  be  solved  by  using  multiple 
methods  to  give  important  information  such  as  the  time  and  distance  a  given  ray  travels 
[19].  Probably  the  most  important  result  of  the  ray  theory  though  is  Snell’s  Law. 


2,  Snell’s  Law 


Snell’s  law  describes  how  a  ray  propagates  in  a  medium  with  changing  sound 
speed.  This  is  especially  helpful  in  modeling  rays  traveling  through  the  stratified  ocean 
since  the  sound  speed  in  water  changes  due  to  pressure,  temperature,  and  salinity.  We 
define  the  sound  speed  of  one  region  to  be  and  the  sound  speed  in  a  second  region  to 

be  Cj .  If  we  know  a  ray’s  angle  of  incidence  onto  the  second  medium  referenced  to 
the  horizontal  then  we  know  the  angle  9^  in  which  the  refracted  ray  will  propagate  going 
through  the  second  region  [18]. 


(5.6) 


cos{9^^  cos{92) 

Cj  Cj 


Snell’s  law  would  seem  to  only  work  between  regions  of  constant  sound  speed. 
However,  any  region  with  non-constant  sound  speed  can  be  considered  as  an  infinite 
amount  of  small  layers  with  “constant”  sound  speeds.  It  can  even  be  shown  that  as  long 
as  the  changing  sound  speed  is  linear  in  the  region  of  interest,  a  ray  will  propagate  in  a 
circular  path  with  radius  R  related  to  the  sound  speed  at  the  location  of  the  ray’s  source, 

Cq  ,  and  the  slope  of  the  sound  speed  profile,  a  [18]. 


(5.7)  R  =  -^ 

a 

Even  with  these  results  it  would  still  be  hard  to  calculate  ray  paths  for  media  in 

which  the  sound  speed  is  not  varying  linearly.  However,  computer  programs  can  use  the 

information  we  just  derived  to  compute  ray  paths.  As  mentioned  before,  the  program 
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numerically  quantizes  the  region  of  interest  into  many  small  layers  of  constant  sound 
speed  in  order  to  approximate  the  ray  path  by  using  Snell’s  Law. 

B,  PROPAGATION  MODELING  USING  BELLHOP 

Bellhop  is  a  ray-tracing  program  that  models  acoustic  pressure  fields  in  any  ocean 
environment.  Given  a  sound  speed  profile,  boundary  conditions,  bathymetry,  and  other 
information  about  the  environment,  it  can  produce  graphs  of  transmission  loss,  eigenrays, 
ray  arrivals,  and  received  time  series  [20]. 

My  goal  of  using  Bellhop  was  to  identify  how  much  of  a  difference  the  depth  of 
the  sensor  would  impact  communications  to  a  sea  floor  node.  The  depth  of  these  sea  floor 
nodes  is  as  much  as  3000  meters.  I  used  the  Munk  deep  sea  sound-speed  profile  to 
represent  a  deep-water  ocean  environment.  The  Munk  profile  is  given  by  the  following 
equation,  where  c(z)  is  the  sound  speed  as  a  function  of  depth  [21]. 


(5.8) 


c(z)  =  1500 


1.0  +  e 


(z-l  +  e“^) 


(5.9) 

(5.10) 


The  scaled  depth  z  and  the  constant  e  are  given  the  following  values. 

6  =  0.00737 

~  2(z-1300) 

1300 

The  resulting  sound  speed  profile  is  shown  in  Figure  13. 
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Munk  Sound  Speed  Profile 


Compressional  sound  speed  (m/s) 

Figure  13.  A  3000  m  Munk  sound-speed  profde. 


Next,  I  used  Bellhop  to  generate  ray  traees  for  a  sound  souree  loeated  at 
0.25  meters  depth  and  then  at  7.0  meters  depth.  Onee  again,  the  goal  was  to  see  if  there 
was  an  appreeiable  amount  of  ray  path  differenee  between  the  two  locations.  The  Wave 
Glider  can  maintain  a  watch  circle  around  a  node  of  about  25  meters,  but  to  be  safe  the 
results  of  the  ray  trace  were  graphed  all  the  way  up  to  1000  meters  to  simulate 
communications  between  the  Wave  Glider  and  a  bottom  node  that  was  farther  away. 
Figure  14  is  a  graph  of  the  ray  tracing  done  with  a  source  depth  of  0.25  meters  and  Figure 
15  has  the  source  depth  at  7  meters.  The  minimum  and  maximum  launch  angles  relative 
to  the  horizontal  axis  were  55  and  89  degrees,  respectively. 
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Bellhop  0.25  meter  Depth  Modem  Simulation  1=1 4000Hz 


Figure  14.  Ray  tracing  for  a  source  at  0.25  m  depth  to  a  range  of  1000  m 


Bellhop  7  meter  Depth  Modem  Simulation  f=14000Hz 


Figure  15.  Ray  tracing  for  a  source  at  7.0  m  depth  to  a  range  of  1000  m 
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As  we  can  see  from  the  graphs,  the  depth  made  very  little  difference  with  respect 
to  the  propagation  paths  of  acoustic  rays  from  the  modem.  Next,  I  checked  if  the  depth 
would  make  a  difference  in  case  the  modem  on  the  Wave  Glider  was  to  communicate 
with  a  node  that  was  much  further  away  than  the  closest  node.  I  extended  the  maximum 
range  calculated  by  Bellhop.  Figure  16  is  a  graph  of  the  ray  tracing  done  with  a  source 
depth  of  0.25  meters  and  Figure  17  with  a  source  depth  of  7  meters.  Once  again,  the 
frequency  used  was  14kHz.  The  minimum  launch  angle  was  -10  degrees  and  the 
maximum  launch  angle  was  10  degrees  relative  to  the  horizontal  axis. 

Bellhop  0.25  meter  Depth  Modem  Simulation  1=14000Hz 
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Bellhop  7  meter  Depth  Modem  Simulation  f=14000Hz 


Although  it  is  hard  to  see,  if  we  compare  the  two  graphs  closely,  we  can  see  that 
when  the  modem  is  at  7  meters  of  depth  it  is  able  to  transmit  a  signal  a  lihle  bit  further 
without  reflecting  off  of  the  bottom.  This  makes  sense  since  a  deeper  source  would  be 
able  to  take  advantage  of  the  downward  refracting  sound  speed  profile  more  than  one 
closer  to  the  surface.  This  means  that  if  we  desire  the  Wave  Glider  modem  to 
communicate  to  further  away  nodes,  it  would  be  advantageous  to  place  the  modem  either 
on  the  glider  portion  or  in  a  tow. 

C.  SHADOWING/INTERFERENCE  OF  GLIDER  BODY 

Some  interference  of  the  wave  glider  itself  could  occur  if  a  piece  of  the  vehicle 
were  to  lie  in  between  the  modem  and  the  sea  floor  node  that  it  is  communicating  with. 
This  problem  could  only  really  exist  if  the  modem  was  placed  on  the  float  portion  of  the 
vehicle  since  the  glider  body  is  the  only  part  of  the  vehicle  that  could  lie  between  the 
modem  and  the  node.  The  glider  portion  of  the  vehicle  tends  to  lead  the  float  portion  by 
0.55  meters  [6].  Using  simple  trigonometry  we  can  calculate  the  horizontal  distance  the 
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node  would  have  to  be  from  the  float  in  order  to  cause  an  issue.  Figure  18  shows  the 
geometry  of  the  float,  glider,  and  sea  floor  node.  It  was  found  that  the  node  would  have  to 
be  235  meters  away  from  the  float. 


Figure  18.  The  distances  involved  with  the  glider  interfering  with  the  communication 

of  a  float  mounted  modem 

At  this  distance  and  depth,  the  propagation  angle  of  the  signal  would  be  85 
degrees  relative  to  the  horizontal.  In  this  sound  speed  profde,  the  signal  would  propagate 
in  a  relatively  straight  line.  Also  if  the  Wave  Glider  wanted  to  communicate  to  further 
away  nodes,  the  distances  involved  would  be  much  greater  than  235  meters.  It  should  be 
noted  however,  that  the  blockage  problem  can  be  completely  avoided  if  the  modem  were 
to  be  placed  on  the  glider  or  on  a  tow  body. 

D.  SUMMARY 

In  all  aspects  of  our  acoustic  considerations,  it  would  be  more  beneficial  to  have 
the  modem  away  from  the  surface.  If  the  Wave  Glider  modem  is  to  communicate  with  far 
away  nodes,  a  glider  or  tow  body  mounted  modem  could  communicate  to  longer 
distances.  Also  a  glider  or  tow-body  mounted  modem  would  completely  avoid 
interference  from  the  glider  portion  of  the  vehicle. 
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VI.  NOISE 


A.  BACKGROUND 

1.  Noise  Overview 

The  amount  of  noise  assoeiated  with  the  eandidate  modem  loeations  is  an 
important  eonsideration.  If  noise  dominates  eommunieation  in  a  partieular  integration 
seheme,  then  the  signal-to-noise  ratio  (SNR)  is  degraded  and  the  aeoustie  modem  will  not 
be  able  to  effeetively  relay  information  without  errors.  The  main  eriteria  that  we  are 
looking  for  in  examining  noise  is  a  high  SNR.  We  ean  use  the  passive  sonar  equation  to 
illustrate  this. 

(6.1)  SL-TL>NL-DI  +  DT 

SL  is  the  souree  level,  TL  is  the  transmission  loss,  NL  is  the  noise  level,  DI  is  the 
direetivity  index,  and  DT  is  the  deteetion  threshold  expressed  in  dB  units.  If  we  rearrange 
this  equation  differently,  we  ean  more  elearly  see  the  SNR. 

(6.2)  SL-TL-(NL-DI)>DT 

Sinee  these  are  dB  quantities,  additions  and  subtraetions  eorrespond  to  multiplieations 
and  divisions,  respeetively.  Here  we  ean  see  that  the  signal  transmitted  to  the  modem 
after  aeeounting  for  transmission  loss  must  be  greater  than  the  noise  by  a  eertain  amount 
referred  to  as  our  deteetion  threshold  [18]. 

Noise  ean  be  separated  into  different  eategories.  An  analysis  of  the  ambient  noise 
and  the  self  noise  assoeiated  with  eaeh  loeation  on  the  vehiele  follows. 

2.  Ambient  Noise 

Ambient  noise  is  defined  as  the  noise  that  is  not  due  to  a  hydrophone  itself  and/or 
its  manner  of  mounting.  In  short,  it  is  the  noise  due  to  the  oeean  environment  itself. 
Ambient  noise  ean  eome  from  different  sourees  in  the  oeean  and  these  different  sourees 
dominate  at  different  frequeney  bands  and  vary  with  weather  and  sea  state  eonditions. 
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Specifically,  we  are  interested  in  noise  in  the  spectral  region  around  lOkHz,  since  the 
modems  operate  at  around  10  kHz.  Figure  19  shows  that  shipping  noise  is  not  of  concern 
above  1  kHz  [18]. 


Figure  19.  Sources  of  ambient  noise  in  the  ocean  and  their  corresponding  frequency 

ranges. 

Most  sources  of  the  noise  in  this  band  originate  at  the  sea  surface  and  experience 
attenuation  with  depth,  although  usually  the  distance  scales  that  people  are  have  studied 
are  much  larger  than  the  7  meters  we  are  looking  to  differentiate.  Results  also  vary 
greatly  with  the  geographic  location  in  which  the  vehicle  is  operating  in. 

Some  sources  of  ambient  noise  are  completely  unavoidable  regardless  of  where 
the  modem  is  placed,  or  vary  only  in  small  amounts  relative  to  the  depth  of  the  modem. 
Examples  include  seismic  disturbances,  oceanic  turbulence,  thermal  noise,  and  ship 
traffic.  We  will  only  examine  the  sources  of  noise  that  vary  with  the  location  of  the 
modem  on  the  Wave  Glider,  since  this  is  the  only  variable  we  have  control  over  [18]. 

B,  ANALYSIS 

1,  Noise  Due  to  Surface  Waves 

Noise  from  surface  waves  can  affect  the  0.5  to  25  kHz  band,  and  is  caused  by  the 

sea  state  and/or  the  wind  forces  on  the  surface  of  the  water.  Better  correlations  are  known 
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for  wind  speed  than  sea  state  simply  beeause  sea  state  is  difficult  to  quantify.  The  exact 
process  by  which  surface  waves  generate  noise  is  not  completely  understood.  However,  a 
theory  has  been  produced  relating  the  period  of  surface  waves  to  how  quickly  their 
pressure  changes  attenuate  with  respect  to  depth.  This  attenuation  can  be  seen  in 
Figure  20  [18]. 


Figure  20.  Attenuation  of  pressure  changes  due  to  surfaces  waves  at  a  given  wave 

period  and  depth 

As  we  can  see  from  the  graph,  the  pressure  due  to  surface  waves  attenuates  quite 
rapidly  as  the  depth  is  increased  and  as  the  period  of  the  waves  is  decreased.  At  the 
frequencies  we  are  interested  in,  the  noise  due  to  surface  waves  will  attenuate  quickly 
with  depth.  Therefore,  it  is  advantageous  to  place  the  acoustic  modem  at  a  depth  of 
7  meters  as  opposed  to  closer  to  the  surface  in  order  to  benefit  from  the  rapid  attenuation 
of  the  surface  wave  noise  with  depth. 
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2. 


Noise  Due  to  Rain 


Rain  can  contribute  an  appreciable  amount  to  the  ambient  noise  speetrum  in  the 
frequeneies  that  we  are  interested  in.  The  magnitude  of  the  effect  can  be  attributed  to  the 
rate  of  the  rainfall  and  the  depth  of  the  modem.  Figure  21  shows  the  level  of  noise  in  the 
frequencies  of  interest  [18]. 


Figure  21 .  Speetrum  levels  of  noise  attributed  to  rain 

If  we  consider  eaeh  drop  of  water  hitting  the  surfaee  of  water  to  be  a 
hemispherieal  souree,  we  ean  get  an  estimate  to  the  transmission  loss  that  the  noise  would 
experienee  going  to  an  aeoustie  modem  loeated  7  meters  below  the  surfaee.  Normally, 
the  surface  area  of  a  spherieally  spreading  source  is  given  by  the  surfaee  area  of  a  sphere 
of  radius  R.  However,  we  are  dealing  with  a  surfaee  area  of  a  hemisphere,  so  the  surface 
area  is  halved. 

(6.3)  -  Inr^ 

The  power  passing  through  this  surfaee  is  given  by  the  integration  of  the  intensity 
of  the  wave  over  this  area. 

(6.4) 

•’  PoC 
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The  total  power  remains  constant  as  the  sound  passes  through  a  surface  r  meters 


away. 

(6.5) 


n  = 


PoC 


PoC 


Therefore,  the  following  ratio  is  true. 


(6.6) 


Prmsi^) 


Transmission  loss  is  defined  as 


(6.7) 


7Z  =  201og 


^  n\ ^ 
Prmsi^) 

^PrmAf), 


201og(r) 


where  the  transmission  loss  is  expressed  in  dB  referenced  to  1  meter.  Coincidentally,  this 
is  the  same  as  if  we  assumed  spherical  spreading.  If  we  plug  in  a  distance  of  7  meters  into 
this  formula,  we  see  that  we  get  a  transmission  loss  of  16.9  dB  re  Im.  Thus,  if  there  is 
rain  in  the  environment  that  the  Wave  Glider  is  operating  in,  the  noise  due  to  the  rain 
would  be  greatly  diminished  if  the  modem  is  placed  either  on  the  lower  body  or  a  tow 
body. 


3,  Self-Noise  of  the  Wave  Glider 

A  full  study  of  the  self-noise  due  to  the  Wave  Glider  itself  has  not  been 
completed.  Due  to  the  propulsion  system  of  the  vehicle  being  entirely  wave  actuated 
though,  the  noise  of  the  vehicle  would  only  be  mechanical  in  nature.  The  only  moving 
parts  would  exist  on  the  lower  body,  and  the  noise  due  to  these  components  would  be 
proportional  to  the  speed  of  the  vehicle.  The  frequency  band  would  also  be  related  to  the 
period  of  the  waves  and  the  sea  state.  Complete  avoidance  of  the  self-noise  of  the 
propulsive  mechanisms  of  the  vehicle  though  is  recommended  and  could  be  achieved  by 
either  placing  the  modem  on  the  float  portion  of  the  vehicle  or  the  tow  body. 
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C.  SUMMARY 

In  all  cases  of  our  noise  analysis,  the  tow  body  location  is  a  favorable  alternative 
to  both  the  float  and  lower  body  loeations.  The  tow  body  would  avoid  more  of  the  noise 
at  the  surface  than  if  it  was  mounted  on  the  float  portion,  and  it  would  avoid  more  of  the 
self-noise  of  the  vehicle  than  if  it  was  mounted  on  the  lower  body.  Once  again,  the 
differences  in  distance  that  we  are  analyzing  are  much  smaller  than  many  of  the  studies  in 
noise  that  have  been  done,  but  regardless  of  how  small  the  distance  is,  any  advantage  that 
the  modem  would  get  in  SNR  should  be  considered. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Although  the  estimated  speed  penalties  seem  low,  the  relationship  between  the 
values  seems  reasonable.  In  terms  of  minimizing  drag,  either  the  float  or  a  tow  body  are 
viable  options.  For  our  aeoustie  analysis,  the  glider  or  tow  body  mounts  would  be  able  to 
communicate  with  further  away  nodes  and  would  avoid  the  problem  of  possible 
interference  with  the  glider  body.  Lastly,  a  float  mounted  modem  would  suffer  from 
possible  surface  noise,  while  a  glider  mounted  modem  would  be  subject  to  the  self-noise 
of  the  propulsive  fins  on  the  glider.  A  modem  mounted  on  the  tow  body  would  minimize 
these  noise  sources. 

From  the  results  in  this  study,  a  tow  body  would  be  the  preferable  location 
because  of  our  acoustic  and  noise  related  considerations.  The  drag  is  worse  than  the  float 
mounted  case;  however,  since  we  are  integrating  an  acoustic  modem  it  is  more  important 
that  the  modem  is  able  to  communicate  effectively.  The  tow  body  mount  also  avoids 
possible  interference  with  the  propulsive  mechanisms  of  a  glider-mounted  modem. 

The  results  in  this  thesis  seem  to  agree  with  other  tests  that  have  been  done  on 
acoustic  modem  integration  with  the  Wave  Glider.  In  one  particular  test,  a  modem’s 
performance  was  compared  in  different  positions  on  the  Wave  Glider.  One  integration 
was  on  the  glider  body  and  the  other  was  on  the  float  body.  In  both  of  these  cases,  the 
modem  demonstrated  effective  communications  with  its  target  to  within  2.5  km. 
However,  it  was  noted  that  the  speed  of  the  vehicle  was  impacted  heavily  when 
integrating  the  modem  with  the  glider  body  [22]. 

B,  RECOMMENDATIONS 

This  analysis  was  all  theory  based  and  many  approximations  were  made  in 
deriving  the  results,  especially  those  of  the  drag  analysis.  Obtaining  more  data  on  the 
relationship  between  Wave  Glider  velocity  and  wave  height  would  make  calculating  the 
power  available  to  the  Wave  Glider  more  accurate.  It  would  also  be  extremely  helpful  for 
this  study  to  acquire  velocity  and  drag  data  for  the  Wave  Glider  from  flow  tank 
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experiments.  Lastly,  quantifying  the  self-noise  of  the  Wave  Glider  in  different  sea  state 
conditions  would  help  identify  the  impact  it  would  have  on  a  receiver  located  on  the 
Glider. 

This  thesis  did  not  examine  the  possibility  of  separating  the  transmitter  and 
receiver  portions  of  the  modem  to  different  locations.  An  in-depth  analysis  of  this 
integration  scheme  should  be  considered  in  future  studies. 
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